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We present a finite-size scaling analysis of the droplet condensation-evaporation transition of a 
lattice gas (in two and three dimensions) and a Lennard-Jones gas (in three dimensions) at fixed 
density. Parallel multicanonical simulations allow sampling of the required system sizes with precise 
equilibrium estimates. In the limit of large systems, we verify the theoretical leading-order scaling 
prediction for both the transition temperature and the finite-size rounding. In addition, we present 
an emerging intermediate scaling regime, consistent in all considered cases and with similar recent 
observations for polymer aggregation. While the Intermediate regime locally may show a different 
effective scaling, we show that it is a gradual crossover to the large-system scaling behavior by 
including empirical higher-order corrections. This implies that care has to be taken when considering 
scaling ranges, possibly leading to completely wrong predictions for the thermodynamic limit. In 
this study, we consider a crossing of the phase boundary orthogonal to the usual fixed temperature 
studies. We show that this is an equivalent approach and, under certain conditions, may show 
smaller finite-size corrections. 


I. INTRODUCTION 

The limit of large systems, or the thermodynamic 
limit, is of general interest when studying finite systems. 
It is this limit, that is commonly accessible to experi¬ 
ments. Recent developments in experimental and sim¬ 
ulation techniques open the door to mesoscopic length 
scales [l|. Here, general concepts of finite-size scaling 
may be tested and verified. Considering, for example, 
the equilibrium properties of several homopolymers, the 
aggregation transition temperature was shown to exhibit 
systematic finite-size effects that deviate from the ex¬ 
pected behavior for particle gas systems [^. It was ar¬ 
gued that the observed scaling is valid for intermediate 
system sizes, where the aggregate includes most of the 
polymers, but in the limit of increasing polymer number 
the particle picture should be recovered. On the other 
hand, it is a reasonable assumption that the intermediate 
regime should be also apparent for particle gas conden¬ 
sation, which will be a main focus of the present paper. 

First-order phase transitions in spin systems are usu¬ 
ally separating homogeneous phases, which leads in gen¬ 
eral to finite-size scaling corrections on the order of the 
inverse system volume. Exceptions occur for example 
if the low-temperature phase is exponentially degener¬ 
ated 3 . The general situation changes when we consider 
a phase transition between a homogeneous phase and a 
mixed phase - as in the case of the droplet condensation- 
evaporation transition separating a supersaturated gas 
phase and a mixed phase consisting of a single droplet 
in equilibrium with surrounding vapor. There exists a 
large amount of theoretical literature on this topic about 
the leading-order scaling behavior of large systems [3- 

with origins already in the 80s Q- Usually, a sys¬ 
tem at fixed temperature is considered, estimating the 
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transition density at which a single macroscopic droplet 
forms. In the limit of very large systems, this transi¬ 
tion density is supposed to have finite-size corrections of 
the order while the rounding of the transition 

should scale with The former was numerically 

verified to leading-order at fixed temperature for a lat¬ 
tice gas [ilM3. Several studies of a three-dimensional 
Lennard-Jones gas pj] - [l^ also verified the finite-size cor¬ 
rections of the transition density. For a fixed density, the 
scaling of the transition temperature was demonstrated 
for the two- and three-dimensional lattice gas 0 , while 
the scaling exponent of the finite-size rounding could not 
be verified. There exist complementary microcanonical 
studies of the Ising model at fixed density (or magneti¬ 
zation) [l^, [T^ , demonstrating that the occurring tran¬ 
sition shows signatures of a first-order transition. 

The present study aims to fill the gap with respect 
to the finite-size scaling of the droplet condensation- 
evaporation transition at fixed density and the gradual 
crossover from an effective intermediate regime to the 
asymptotic large-system limit. In order to draw gen¬ 
eral conclusions, we consider a lattice-gas model in two 
and three dimensions as well as a Lennard-Jones gas in 
three dimensions. We will show that the proposed scal¬ 
ing behavior of the transition temperature and round¬ 
ing is numerically recaptured in the limit of large (but 
achievable) system sizes for all dimensions and models 
considered. Where possible, we will directly compare to 
analytic solutions or low-temperature series expansions. 
In addition, we discuss an emerging intermediate scaling 
regime, best visible in the rounding of the transition and 
consistent with recent observations of finite-size effects 
in polymer aggregation [2|. This regime with effective 
local scaling behavior and the gradual crossover to the 
large-system regime may be described by the finite-size 
scaling behavior when including empirical higher-order 
correction terms. 

The remaining part of this paper is organized as fol¬ 
lows. In Sec. m we briefly recapture the main leading- 
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FIG. 1. Sketch of the infinite system-size transition (solid 
line) together with the finite-size scaling directions in either 
the density {T fixed) or temperature (p fixed). The transition 
line may be understood as To{p) or similarly po{T). At the 
crossing point of both schemes a finite system of size V may 
be constructed for which (p, T) = {pc ,Tc)v corresponds to 
the finite-size condensation-evaporation transition. 

order results from the literature at fixed temperature 
and convert them to the scenario of fixed density. Af¬ 
ter describing our models and methods in Sec. Illli we 
will present our results in Sec. IIVI and finish with the 
conclusions in Sec. 0 

II. THEORY 

We consider a d-dimensional liquid-vapor system of 
N particles in a (periodic) box of volume V, where for 
sufficiently dilute systems the condensation-evaporation 
transition separates a homogeneous supersaturated gas 
phase from a mixed phase of a droplet in equilibrium 
with surrounding bulk gas. At fixed density, this occurs 
at the condensation-evaporation transition temperature 
Tc- However, the usual description of this problem is 
in the language of a grand-canonical scheme, considering 
a fixed temperature and variable density (or in case of 
a lattice gas equivalently the magnetization of the Ising 
model) with the finite-size transition density pc- This is 
orthogonal to the scenario of a fixed density with variable 
temperature. Figure [T] shows a sketch of the infinite-size 
phase boundary Tq{p) or po{T) and the different scenar¬ 
ios crossing it for finite system sizes. It needs to be men¬ 
tioned that both schemes are working in the canonical 
ensemble for each point (p, T). So in fact, any finite-size 
transition point {pc, Tc)v belongs to one fixed-T and one 
fixed-p scheme, simultaneously. The same holds for any 
canonical function /(p, T). Thus, the orthogonal crossing 
schemes are equivalent and we may translate a functional 
dependence f{p,T)V°‘ = 1 from one scheme to another 
by a Taylor series expansion. Then, expanding around 
some T* yields 

= /(p, T*) + f{p, T*){T - T*) + ... , (1) 

which may be solved for T. The remaining task is to 
identify and evaluate suitable functional dependencies. 


Using this procedure, we will recapture for a fixed 
density p = N/V the finite-size correction to the tran¬ 
sition temperature Tc — Tq oc and the scaling 

of the finite-size rounding (or width of the finite-system 
transition region) AT = T — (x This is 

directly related to the linear extension of the droplet 
R oc which becomes the relevant length scale 

for large system sizes. 


A. Finite-size scaling of the droplet 
condensation-evaporation transition 

Biskup et al. Q showed for a d-dimensional liquid- 
vapor system in equilibrium a vanishing probability of 
intermediate-sized droplets. This was proven rigorously 
for the lattice-gas interpretation of the 2D Ising model 
and justifies the restriction of the discussion to two rele¬ 
vant phases: the homogeneous gas phase and the mixed 
phase of a single macroscopic droplet in equilibrium with 
a bulk gas surrounding it. At coexistence both phases are 
equally probable and in general (finite volume) the dom¬ 
inance of either phase is the solution of a saddle point 
problem (see also Ref. Q). Choosing a fixed tempera¬ 
ture T, we consider a system at supersaturated density 
p = N/V with particle excess SN = N — Nq = SpV, 
where po = Nq/V is the equilibrium density of the 
infinite-system gas phase at this temperature. In fact, 
the particle excess is temperature dependent, see Fig. [T] 
A fraction A € [0,1] of this particle excess is considered 
to be in the single macroscopic droplet, i.e., JYd = X6N, 
while the remaining excess is in the gas phase such that 
SN = i5Yd + SNq. Up to higher-order corrections from 
surface and density fluctuations as well as translational 
entropy of the droplet motion, Biskup et al. Q intro¬ 
duced a rescaled (size-independent) density 

, (pL-po)^(^iV)^ 

2krw V ’ ^ 

with the infinite-system liquid and gas densities pn and 
Po, the Wulff-shape surface free-energy per unit volume 
Tw and the (reduced) isothermal compressibility k, see 
also Ref. 0- Here, only SN is variable and the remain¬ 
ing parameters are constants encoding all system-specific 
details, like the droplet shape. In the limit of large sys¬ 
tem sizes, they calculate the fraction of particles in the 
largest droplet A(A) as a “universal” function of A. At 

the transition density A^ = ^ system 

changes between the vapor phase (A = 0) and the mixed 
phase (Ac = 2/(d -I-1)). 

Now consider a fixed density p with varying tem¬ 
perature T. Then the system-specific constants be¬ 
come functions of the temperature, namely Pi{T), k{T), 
tw(T) and SN{T) = {p — po{T))V. The temperature- 
dependent particle excess may be understood by the fact 
that for decreasing temperature (fixed density) the dif¬ 
ference Sp increases and with it the supersaturation, see 
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also Fig. [TJ We may rewrite Eq. m as 


A^V-^ = f{p,T), (3) 


where we identify 


f{p,T) 


p-po{T) / (pL(T)-po(r))^ \^ 
pUT)-Po{T)\ 2/t(r)rw(T) J 


At the condensation-evaporation transition, A = Ac is 
constant and the left-hand side of Eq. m is depend¬ 
ing only on the system size V. Then, for a fixed fi¬ 
nite system size, a suitable combination of T and p 
solves Eq. (|31) yielding the finite-size transition point at 
{p,T) = (pc,Tc)v in Fig. [TJ This transition point may 
be obtained either numerically exact or by a Taylor ex¬ 
pansion. Keeping p = NjV constant, we proceed by 
expanding f{p,T) in Eq. (JT]) around the infinite-system 
transition temperature Tg, where po{To) = p and thus 
f{p,To) = 0. Solving for the finite-size transition tem¬ 
perature T = Tc yields 


Tr = Tn 


A 




np,To) 


V <i+i -I- O (F 'i+i 


(v 


( 5 ) 


In terms of the number of particles this means to first 
order 


Tc-To(x N~^. ( 6 ) 

We need to emphasize that Eq. ((Sj) is only the leading- 
order result. Numeric tests at fixed temperature already 
showed apparent higher-order corrections 1M3- 

For the Ising model, we can go a little further and 
evaluate the leading finite-size scaling explicitly. To this 
end, we relate the magnetization and the density via 
TO = 1 — 2p and identify the spontaneous magnetization 
Too = 1 — 2po as well as, from the symmetry of the Ising 
model, pl = 1 ~ Po- Iii addition, we find the magnetic 
susceptibility x = k and = 4rw, due to the shift in 
the energy scale when exploiting the equivalence of the 
Ising model and the lattice gas model 0- Thus, we 
obtain 




1 / TO W 2too(TI")2 \ <^+1 

2 


( 7 ) 

with = 4T when rewriting the Hamiltonian according 
to the definition in Sec. nni and Ref. 0- 

In two dimensions, the involved quantities are known 
analytically or up to arbitrary precision: toq is de¬ 
scribed by the Onsager-Yang equation [13, X is ob¬ 
tained from sufficiently long series expansions [Tsl - I^ 
and = 2\/W can be obtained from the volume of 
the Wulff plot W [2lj. For a collection of equations we 
refer to Ref. Q. We numerically evaluated Eq. (Q in 
Eq. ((Sj) , fixing the density and volume and solving with a 
bisection algorithm for the corresponding transition tem¬ 
perature. We will compare this result to the numerical 


finite-size results later in Fig. |6| (top) and refer to it as 
the full solution of Eq. ©. 

In three dimensions, we may make use of low- 
temperature series expansions of the spontaneous mag¬ 
netization [ 23 . This allows us to estimate the 
infinite-system transition temperature Tq by solving 
p — po(T) = 0. We will compare these results to the 
finite-size scaling fits of our numerical data and refer to 
it as the solution from low-temperature series expansion 
in Fig. ini (center). 

The analytic results may be used to recapture the 
scaling of the droplet size at condensation, where the 
fraction Ac will be in the largest droplet. For the vol¬ 
ume of the droplet at condensation it follows 5Vu = 
(pl — po)~^\c5Nc [ 1 |, where the total particle excess at 
condensation 5Nc encodes the shape of the droplet (see 
Eq. (l2|)). In general, the volume of an ideal droplet may 
be expressed by 5Vd = SdR'^, where Sd is a geomet¬ 
ric shape factor that allows to describe both spherical 

and cubic droplets, where the latter may occur in lattice 
systems below the roughening transition M- Equating 
both droplet volumes, inserting SNc from Eq. ([2|), and 
solving for the radius yields 

R= {Sj^SVoy (8) 

consistent with results in the literature H, 01 ■ For fixed 
density, this is equivalent to i? oc Moreover, 

the leading finite-size scaling corrections in Eq. may 
be expressed in powers of R~^, i.e., Tc — Tqoc. R~^. 


B. Rounding of the droplet 
condensation-evaporation transition 


In analogy, one may argue that the rounding of the 
transition at fixed density should scale with the system 
size in the same way as the rounding of the transition 
density (or magnetization) at fixed temperature. The lat¬ 
ter was derived by Binder [3 using a phenomenological 
theory. He predicted a finite-size rounding proportional 
to Starting from a two-state approximation 

the discussion may be reduced to a homogeneous and an 
inhomogeneous phase. Then, the rounding of the tran¬ 
sition may be related to a characteristic width fiAF of 
order unity, where AF is the free-energy difference be¬ 
tween both phases. 

Expanding the free-energy difference /3AF around 
the finite-size transition temperature Tc, recalling that 
{d/dp)PF = T, yields 


/3AT = {f3AF)\^^ 



(T-Tc) + ... . (9) 


The free-energy difference vanishes at Tc in the limit 
of large system sizes, considering that both phases con¬ 
tribute with equal probability. In the gas phase, the par¬ 
ticles may be considered non-interacting. Thus, the en¬ 
ergy difference is dominated by the droplet energy, which 









4 



FIG. 2. Sketch of the particle gas models under investiga¬ 
tion: The discrete lattice gas model (left) in two and three 
dimensions and the continuous Lennard-Jones gas model with 
domain decomposition (right) in three dimensions. 


by only one particle. Introducing a state-variable for each 
site, rii £ {0,1}, the Hamiltonian may be written as 

— ( 11 ) 

equivalent to the Ising model (si £ {if}) at fixed mag¬ 
netization if TT-i = i {si -I-1). Here, the density p = N/V 
can only be adjusted approximately by occupying an in¬ 
teger number of sites N = pL‘^ . 


B. Lennard-Jones gas 


depends on the droplet volume (x R‘^. At Tc, Eq. (IS|) 
relates the droplet radius R to the system volume and 
the energy difference may be approximated to leading 
order as AE ^ The finite-size corrections from 

the transition temperature appear merely as corrections 
to the energy difference, such that Eq. (I5|) simplifies to 
f3AF ~ (H'^A'i+i) jk-QT^) AT in the limit of large system 
sizes. The condition |/?AE| ~ 1 yields to leading order 
the rounding width AT oc and in terms of the 

particle number 

AT oc A~^+T. (10) 

Notice that the radius of the droplet R as relevant length 
scale was used in the argumentation such that for large 
system sizes the rounding may be identified as AT = 
T - Tc oc 


III. MODELS AND METHODS 

In order to investigate the universal aspects of the 
condensation-evaporation transition, we employ two dif¬ 
ferent particle gas models, for a sketch see Fig. [2] Bridg¬ 
ing the gap to analytic solutions Q , we consider a lattice 
gas model equivalent to the Ising model at fixed magne¬ 
tization d, [T^. We add to our investigation a Lennard- 
Jones gas model, which — in contrast to the lattice gas — 
is not symmetric with respect to particle-hole exchange 
(see also The boundary conditions are peri¬ 

odic, which is common for a finite-size scaling analysis. 
That way, non-intended interactions with the boundaries 
are avoided. For both models, we consider a finite in¬ 
teraction range smaller than the linear extension of the 
system, which safely allows the use of periodic boundary 
conditions. 


A. Lattice gas 

In the discrete lattice gas (DLG) model, particles are 
described by occupied sites on a two-dimensional (2D) 
square or three-dimensional (3D) cubic lattice of size 
V = L‘^. The interaction is purely short-range, with 
attraction between nearest neighbors and hard-core re¬ 
pulsion by the condition that each site can be occupied 


The Lennard-Jones gas (LJG) model is a continuous 
particle gas model, where in principle all particles inter¬ 
act with each other via the potential 


bij(Dj) = 4e 




( 12 ) 


with e = 1 and a = 2“^/® such that the potential min¬ 
imum is at Tmin = 1- Being consistent with literature 
and reducing the computational demand, we introduce 
a cutoff radius Tc = 2.5a above which the potential is 
zero. The remaining potential is then shifted by lAj (re) 
in order to be continuous: 


KjW = 


4ij(r) - VLj(rc) r<rc 
0 else 


(13) 


That way, we may use a domain decomposition, only 
calculating interactions with particles in neighboring do¬ 
mains. This reduces the computational demand espe¬ 
cially in the gas phase. The Hamiltonian is given by 

= (14) 

i¥^3 

Here, we consider only the three-dimensional case of par¬ 
ticles in a cubic box of size V = with periodic bound¬ 
ary conditions and L = [N/p^/^. 

For the present study, we considered the Lennard- 
Jones gas with the same density as for the lattice gas. 
The Lennard-Jones gas may be reasonably applied to 
non-polar gases, for example Argon (Ar). For this case, 
molecular dynamics simulations were matched with ex¬ 
periments already in the 60s [1^ . Gorresponding param¬ 
eters are a ~ 3.4A and e/fce ~ 120Ar, which serve well for 
an order-of-magnitude comparison. This would lead to 
a real temperature = Te/Zcb. In order to compare 
to the literature boiling point, the density would have 
to be adjusted accordingly, e.g., for Argon, Tboiimg ~ 
87.3K with a gas density pboiimg ~ 5.772g/l at at¬ 

mospheric pressure. The unit length a in our system is 
related to the parametrized length scale as a = 2^/^a. 
Assuming a molar weight of 39.95 x 1.6605 x 
for Argon, this yields the conversion for the density 
^reai ~ p x 1.226 X IQ^g/l. In dimensionless units, the 
experimental boiling (evaporation) temperature and gas 
density yield Tq ~ 0.728 and p « 0.005, respectively. 
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FIG. 3. (Color online) Example of equilibrium droplets at coexistence from Metropolis simulations for the 3D lattice gas with 
N = 10 000 (left) and the 3D Lennard-Jones gas with N — 2 000 (right). 


C. Parallel multicanonical simulations 

The droplet condensation-evaporation transition 
shows a phase coexistence between a supersaturated 
gas phase and a mixed phase consisting of a droplet 
and a remaining gas phase. The resulting free-energy 
barrier i m decreases the probability to pass from 
one phase to the other in the canonical ensemble which 
leads to effects like hysteresis and other effects charac¬ 
teristic for first-order like phase transitions. In order 
to overcome these barriers and to sample the transition 
point with hig h accuracy, we apply multicanonical 
simulations [^, This is a natural choice for the 

fixed-density scheme, while for fixed-temperature a sim¬ 
ilar approach is not straight forward. The Boltzmann 
weight exp(—/Jil), where (3 = (/cbT)“^ with /cb = 1, 
is replaced by an a priori unknown weight function 
which is iteratively modified in order to yield a flat 
histogram This allows the sampling of suppressed 
states between the coexisting phases and produces 
accurate estimates of the observables when reweighting 
around the transition temperature. As a result, we 
obtain a full temperature range from a single simulation 
applying standard time-series and histogram reweighting 
techniques. 

However, this requires to sample a broad energy range 
with sufficiently many tunnel events across the transi¬ 
tion point from high energy to low energy and vice versa. 
We make use of a parallel implementation of the multi¬ 
canonical method |28j . which allows a reduction of com¬ 
putational time. In the cases of lattice gas condensa¬ 
tion, it was tested and shown to speed up the simulation 
by the number of cores used 0- We considered lattice 
sizes up to 1000^ for the 2D and 100^ for the 3D lattice 
gas as well as Lennard-Jones particle systems with up to 
512 particles with a fixed density p = 10“^. The Monte 
Carlo updates in the lattice systems include local particle 
shifts to a nearest-neighbor site combined with particle 


displacements to a random new site. For the Lennard- 
Jones gas, we restrict ourselves to local random particle 
displacements. 


D. Metropolis simulations 

In addition, we applied standard Metropolis simula¬ 
tions in order to sample the droplet phase at coexistence. 
We can use the hnite-size scaling corrections from the fit 
to the transition temperatures obtained with the paral¬ 
lel multicanonical simulations. Then, we make use of 
the first-order nature of the transition causing a suppres¬ 
sion of the intermediate states between the vapor and the 
mixed phase. Preparing a system in the droplet phase, 
the Metropolis simulation thermalizes and only samples 
the droplet phase for a finite time. With increasing sys¬ 
tem size, this finite time gets larger and we may safely 
sample the size of the critical droplet with a lot less statis¬ 
tics than would be required for the multicanonical sim¬ 
ulation. However, we need the multicanonical data of 
suitable system sizes in order to obtain valid estimates of 
the transition temperature. We do not use the Metropo¬ 
lis data to draw conclusions about the transition tem¬ 
perature or the finite-size rounding but instead directly 
measure the critical droplet size. This allows to reach an 
additional order of magnitude in N. An example for the 
3D cases is shown in Fig. [31 


E. Observables 

The focus of this study lies on two observables and 
their thermal derivatives, namely the energy E and the 
fraction of particles in the largest cluster -q = Nd/N. 
The latter is determined by identifying the largest clus¬ 
ter as the maximal number of connected particles Nd . In 
case of the lattice gas, the connected particles are simply 
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FIG. 4. (Color online) Canonical estimates of the fraction 77 of particles in the largest droplet for a 3D Lennard-Jones gas at 
fixed density p = 10~^ (left). Corresponding probability distribution of 77 at Tc for N = 512 (right). 


nearest neighbors. For the Lennard-Jones gas particles 
are defined as connected if < 2a. Thus, 77 is a mea¬ 
sure of the “mass” of the largest droplet and allows one to 
study the homogeneous nucleation transition efficiently. 
An example for the 3D Lennard-Jones gas is given in 
Fig. |4] (left). The consideration of these types of geo¬ 
metric clusters is a safe choice for dilute systems. Note 
that for rather dense lattice systems also the stochas¬ 
tic Swendsen-Wang cluster definition has been consid¬ 
ered [ 2 ^. 

Anticipating a first-order phase transition, we expect 
in the thermodynamic limit a discontinuity in the energy 
as well as in 77 . Hence, the thermal derivatives of the ob¬ 
servables 7^0 = {{EO) — {E){0)) will show a pro¬ 

nounced peak at the transition temperature Tc for finite 
systems, see Fig. [SJ The determination of the peak may 
be done very precisely due to the full temperature range 
from multicanonical simulations. The error is estimated 
by using jackknife error analysis [s^, where we make use 
of the independent parallel production runs combining 
all but one time series for each jackknife bin. In the case 
of the energy, the thermal derivative is related to the 
specific heat Cy = — {E)^) /N. 

Moreover, we will use the mass-fraction 77 in order to 
directly verify the predictions about the scaling of the 
critical droplet size that was assumed and used in the 
theory. We want to focus on the droplet phase and esti¬ 
mate the average fraction of particles in the largest clus¬ 
ter at the transition temperature {f))Tc- To this end, 
we reweight our multicanonical data to droplet size dis¬ 
tributions P{vi) by time-series reweighting, adding each 
reweighting factor to the histogram bin corresponding to 
the observable. An example is shown in Fig. 2] (right). 
Assuming that both phases are at equal weight, we iden¬ 
tify 77min such that Z = j^ dipPirf) = 0.5. Then 

{v)t, = ^[ dr] r]P{r))\j.^ (15) 

is a robust estimator for sufficiently large systems. In 
fact, we directly deal with the number of particles in a 


droplet such that the integrals become discrete sums. Er¬ 
rors are again estimated using the jackknife error anal¬ 
ysis. Systematic errors may arise from imperfect esti¬ 
mates of the coexistence temperature, the definition of 
the droplet boundary and, especially for small systems, if 
the distribution does not show a sufficiently pronounced 
dip between the two phases. 


IV. RESULTS 


The results are presented for a fixed particle density 
p = N/V = 10“^, up to one exemplary case when com¬ 
paring the Lennard-Jones transition temperature to the 
Argon boiling point. This choice of a dilute system en¬ 
sures that the observed transition is really that associated 
with the formation of a droplet and not of a cylinder or 
slap [T^. We have to emphasize in the beginning that 
the lattice model obviously shows discretization effects 
for small systems because this explicit density may not 
be realized for any system size and in general remains 
approximate. We will see, however, that this does not 
influence our main conclusions. 

We consider only a single density, because we want 
to test the leading-order finite-size corrections and scal¬ 
ing exponents. The density will, however, influence the 
finite-size scaling limit. In order to illustrate this, con¬ 
sider the low-temperature series expansion of the 3D 
Ising model. Here, the spontaneous magnetization is 
given to first order as mo = 1 — -|-... . For T —>■ 0 

or po —>■ 0 , this yields for the density po = ^(tuq — 1 ) = 
(1 -|- ...) and to first order for the inverse temper¬ 
ature (3 ~ i In Po • This result may be similarly de¬ 
duced from microcanonical arguments in continuous sys¬ 
tems (see Ref. i)- 
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FIG. 5. (Color online) Canonical estimates of the specific heat 
from parallel multicanonical simulations of the 3D Lennard- 
Jones gas at fixed density p = 10“^. 

A. Finite-size scaling of the droplet 

condensation-evaporation transition temperature tn 

The transition temperature, where the transition- 
related observables show a discontinuity, may be esti¬ 
mated by the peak locations of their thermal derivatives, 
for an example see Fig. [5j In the limit of large systems, 
this transition temperature is predicted to scale as 

Tc = To + o,N 'i+i -|- ..., ( 16 ) 

see Sec. Ill A1 Figure [6] shows numerical results for the 
three considered situations, the lattice gas in two and 
three dimensions and the Lennard-Jones gas in three di¬ 
mensions. The estimated transition temperatures are ob¬ 
tained from the maxima of the specific heat (red pluses) ^ 
as well as the thermal derivative of the fraction of parti¬ 
cles in the largest cluster ^ (green crosses). Both esti¬ 
mates are remarkably similar as expected for first-order 
phase transitions, and hence we only show local fits to 
the estimated transition temperature derived from the 
specific heat. 

In all three cases, we observe for large system sizes a 
proper finite-size scaling behavior according to the pre¬ 
dicted scaling, shown by a good quality of the least- 
squares fits (a reduced per degree of freedom of 
about 1). In accordance with the literature, we ver¬ 
ify that this requires rather large system sizes for the 
leading order fit (dashed dark blue fit). The interme¬ 
diate sized systems may be included in an empirical 
fit with the next higher orders (dotted light blue fit), 
i.e., Tc = To -I- -|- In addition, we 

want to mention that for intermediate system sizes an ef¬ 
fective different scaling behavior may be observed. While 
not completely apparent, it may be justified by suffi¬ 
ciently good x^i that the transition temperature can be 
locally describes by a behavior (dash-dotted or¬ 

ange fit). For details of the individual fits see Table U and 
the following discussion. The scaling of the intermediate 
regime is consistent with studies of flexible homopolymer 



FIG. 6. (Color online) Finite-size scaling of the droplet 
condensation-evaporation transition temperature for the 2D 
lattice, 3D lattice and 3D Lennard-Jones gas from top to bot¬ 
tom. The lattice gas is compared to theoretical predictions 
(see text and also Sec. IIIA|. 

aggregation Q, where a large fraction of the system is 
involved in the formation of the droplet/aggregate. If al¬ 
most all constituents are involved in the transition, then 
the linear extension of the homogeneous, isotropic con¬ 
densate is just - which justifies the assumed scaling 
behavior. This leads to two possible conclusions: Either 
the local scaling function is a polymer property or it is a 
generic property for a small number of polymers, where 
the latter seems to be in agreement with the presented 







































TABLE I. Results of different fit functions to the condensation-evaporation transition temperature with the statistical errors of 
the fits. If no upper range is provided, it refers to A^max = 10 000 for the lattice gas cases (2D and 3D) and A^max = 512 for the 
Lennard-Jones gas. The reference infinite-size value Tq is obtained for the lattice systems from the Onsager solution (2D) Ea 
and from low-temperature series expansions (3D) f^. 


Model 

To 

Range 

■f 

To 


To -t aN- 

Range 

■i/C+i) +bN- 

To 

-2/(d+l) 

To -t 

Range To 


Ref 

To 

2D DLG 

[2500: 

] 

0.39884(3) 

0.4 

[ 400: 

] 

0.3982(1) 

1.8 

[324:900] 

0.3903(2) 

8.0 

0.39882... 

3D DLG 

[2160: 

] 

0.62341(4) 

0.6 

[1663: 

] 

0.6229(4) 

1.9 

[ 68:243] 

0.5840(3) 

1.6 

R 0.622 

3D LJG 

[ 160: 

] 

0.7106(4) 

0.8 

[ 10: 

] 

0.7011(4) 

1.1 

[ 12: 48] 

0.6597(4) 

1.0 



results. In order to explore this observation further, we 
will investigate the size of the critical droplet and the 
rounding of the transition (which in fact provides suit¬ 
able fitting ranges for the transition temperature) in the 
following subsections. 

In the case of the 2D lattice gas, we may compare 
directly to the analytic solution of the infinite system 
and to the full solution of Eq. ([3]) for finite systems de¬ 
rived in Sec. Ill Al The full solution shows large devia¬ 
tions for small systems as expected. With increasing 
system size, however, it starts to describe the finite- 
size scaling approximately. This already gives a hint 
to the choice of a proper htting-range with a leading- 
order scaling behavior. A least-square fit of the leading- 
order for the largest system sizes N > 2500 yields an 
adequate ~ 0-4 and an infinite-size transition tem¬ 
perature To = 0.39884(3). This is consistent with the 
analytic result inverting Onsager’s solution for the mag¬ 
netization [^, with T = r^®/4, 


mo(fo) = 0.98 


l-sinh“^(l/2fo) 



(17) 


The resulting Tq = 0.39882 is shown in Fig. [ 6 ] (top) by 
the arrow. Including the next order and fitting N > 400 
yields a ~ 1.8 with To = 0.3982(1), which slightly de¬ 
viates from the exact solution. This may be taken as a 
hint that our next-order term is only an effective correc¬ 
tion and additional corrections of the same order may be 
apparent. Interesting to compare is also the amplitude 
of the leading-order correction The linear and 

higher order fits yield a = —0.234(1) and a = —0.214(2) 
respectively. This may be compared to the power-series 
expansion of the full solution (see Eq. ([5]) and Eq. ©)■ 
Making use of the analytic solution for mo, the series 
expansion for x up to 300*^ order and the integral so¬ 
lution of Tw (for a list of equations see Ref. [9|), we 
may numerically differentiate Eq. 0 and calculate in 
lattice gas units a = Ac^^p^/^/4/'(to, Tq®) Ri —0.239 in 
decent agreement with the leading-order fit. The previ¬ 
ously mentioned intermediate scaling regime is not very 
prominent for the 2D lattice gas where a least-square fit 
to in the (already small) range N = [324, 900] still 

yields a « 8 . Moreover, the infinite-size extrapolation 
is obviously wrong. It is worth noting, that the inter¬ 
mediate scaling regime in polymer aggregation was 


observed in three dimensions, which suggests that the 
prominence of this regime may depend on the dimension. 

For the 3D lattice gas, we may compare the finite-size 
scaling results of the infinite-size transition temperature 
to low-temperature expansions, see Sec. Ill Al A fit of the 
leading-order scaling behavior to sufficiently large sys¬ 
tems in the range N > 2160 yields Tq = 0.62341(4) with 
Ri 0.6 which is in the vicinity of the (not exact) low- 
temperature expansion Tq r; 0.622 again shown by 
the arrow. Increasing the fit range still yields reason¬ 
ably good fits with higher x^ of the same quality as the 
fit including the next order see also Table H] No¬ 

tice that in the case of the 3D lattice gas also the small 
system sizes seem to coincide with the leading order fit. 
However, considering only an intermediate regime allows 
to fit the behavior with qualitatively good local 

agreement. If the larger system sizes were not present, 
this could be interpreted as the leading order scaling be¬ 
havior. Comparing the fit to this effective ansatz with the 
low-temperature series expansion shows, however, strong 
deviations. 

In the case of the 3D Lennard-Jones gas one may best 
see the arising peculiarities. The leading-order fit for 
N > 160 to yields Tq = 0.7106(4) with r; 0.8 

but shows a clear deviation for small system sizes. In¬ 
cluding the next order for N > 10 yields Tq = 0.7011(4) 
with x^ ~ 1 - 1 , in rough agreement with the leading- 
order result, and recaptures the deviation of the small 
system sizes. This is consistent with results for the same 
Lennard-Jones model PHI. However, considering an 
intermediate regime N = [12 : 48] with the ansatz 
yields a qualitatively good fit with Tq = 0.6597(4) and 
X^ R! 1, which deviates strongly from the fit. 

Again, if the largest system sizes were not present this 
could be interpreted as the finite-size scaling corrections, 
especially if no reference temperature is available. Lo¬ 
cally it seems that this is an intermediate scaling regime 
which is, however, already covered by the theoretically 
predicted scaling behavior including empirical higher- 
order corrections. This shows the necessity to control 
that the available data is really in the expected leading- 
order scaling regime. 

Comparing to the boiling temperature of Argon, we 
have to exemplary consider a density p = 5 x 10“^. Then, 
the leading-order fit yields Tq = 0.6525(3) (not shown 
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FIG. 7. (Color online) Example of the transition rounding 
for the specific heat of the 3D Lennard-Jones gas. The x 
and y axes are rescaled according to the leading-order scaling 
behavior. 

here) which differs from the experimental result Tq 
0.728 (see Sec. IIIIBIl by ^ 10%. While the order-of- 
magnitude is comparable, the difference is expected for 
the truncated Lennard-Jones potential [H. 

The largest systems considered for the Lennard-Jones 
gas included 512 particles in a box of length L Ri 37.1 ~ 
41.7(7. This is a lot smaller than the system sizes con¬ 
sidered in Ref. [l^ (L < lOOcr) at fixed temperature 
T R! 0.68 (for their parameterization) with typical parti¬ 
cle numbers N r 15 800. Still, in this fixed-T approach 
they did not see the predicted scaling behavior of the 
transition density shift but needed to extrapolate an ef¬ 
fective exponent (smaller than —0.89) in order to recover 
the theoretical prediction Similarly, a direct fit 

of the leading-order power-low exponent to the present 
data yields oc already for smaller sys¬ 

tem sizes N > 160. This implies that an orthogonal 
phase boundary crossing may lead in certain situations 
to reduced finite-size corrections and serves as a useful, 
complementary approach. 

B. Finite-size rounding of the droplet 
condensation-evaporation transition 

In the limit of large systems, the rounding of the tran¬ 
sition is predicted to scale as AT oc for fixed 

density, see Sec. Ill Bl Figure [7] shows the rescaled specific 
heat from Fig.[5]for the 3D Lennard-Jones gas. At a first- 
order transition, the maximum of the specific heat should 
scale inverse proportional to the rounding of the transi¬ 
tion . The rounding of the transition is here estimated 
as the half-width of the specific-heat peak, defined as the 
width where Cy > Errors are obtained by jack¬ 

knife error analysis. 

Figure [5] shows the finite-size scaling of the round¬ 
ing for all three considered systems. In all cases, we 
could clearly identify two different scaling regimes: For 
intermediate system sizes one may observe AT oc N~^ 



10 100 1000 10000 


N 



FIG. 8. (Golor online) Finite-size rounding of the droplet 
condensation-evaporation transition for the 2D lattice, 3D lat¬ 
tice and 3D Lennard-Jones gas from top to bottom. The re¬ 
sults are presented for p = 0.01 but are consistent for different 
densities. 

(dashed line) and for large system sizes AT oc 
(solid line). The latter is the expected scaling of the 
rounding as predicted by theory i[ll and reformulated 
in Sec. Ill Bl Thus, having only data available in the in¬ 
termediate regime would suggest a wrong hnite-size scal¬ 
ing behavior consistent with the intermediate regime for 
the transition temperature. However, the intermediate 
regime is nicely htted by considering the predicted scaling 
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behavior for large system sizes and empirically including 
the next higher-order corrections (dash-dotted light blue 
fit), i.e., AT = -h b'. 

In the case of the 2D lattice gas, smallest systems show 
no systematic behavior which explains the large devia¬ 
tions in Fig. |6l For the intermediate regime one can see 
a direct particle or volume dependence (p fixed). The 
onset of the large-system regime is consistent with the 
scaling range in the transition temperature that approx¬ 
imately coincides with the full solution (Fig. [B]). A fit 
including higher order corrections yields a Ri 0.7 in¬ 
cluding already system sizes N > 324 and thus including 
the intermediate regime. 

For the 3D lattice gas and the Lennard-Jones gas, the 
intermediate regime is apparent for quite small systems 
and also shows the A-proportionality up to the crossover 
to the large-system regime. Again, the crossover is 
consistent with a good choice of a leading-order fitting 
range for the finite-size scaling of the transition tempera¬ 
ture. The fit to the rounding of the transition, including 
higher-order corrections, allows to include system sizes 
N > 175 with Ki 2.3 for the lattice gas and A > 16 
with x^ ~ 0.9 for the Lennard-Jones gas case. Again, 
this includes (parts of) the intermediate regime by con¬ 
sidering empirical higher-order corrections. 

Previous studies of the lattice gas in two and three 
dimensions at fixed-p showed significant deviations from 
the predicted exponents for the transition rounding 
using average densities of states from Wang-Landau sim¬ 
ulations. For the 3D case they found an effective scaling 
of the rounding with L~2.45(2) ^ ^ present 

results on the other hand clearly confirm the “large”- 
system scaling behavior. Direct fits of the power-law be¬ 
havior to the largest system sizes of the 3D lattice gas 
model yield effective exponents and 

for N > 1663 and N > 5120 respectively, close to the 
predicted scaling 

The finite-size rounding shows to be a good observable 
to identify the previously noticed intermediate regime. 
The width of the transition is associated to the fluctu¬ 
ations of the system [s^, which should depend on the 
inverse volume of the relevant system size. The relevant 
system size remains the droplet, which scales for large 
systems as R‘^ oc Qn the other hand, for small 

systems the droplet includes a large fraction of the sys¬ 
tem, see the following discussion. Thus the local scaling 
of the intermediate regime is not surprising. Knowing 
the scaling behavior for large systems, however, allows to 
anticipate this gradual crossover with empirical higher- 
order corrections. 


C. Scaling regimes of the critical droplet size 

We identified the leading-order finite-size scaling cor¬ 
rections as powers of the critical droplet radius. This 
requires that the mass of the critical droplet scales in 
the limit of large systems as yve want to test 

this relation by considering the fraction of particles in 




N 



FIG. 9. (Color online) Scaling of the particle fraction in the 
largest droplet p for 2D lattice, 3D lattice and 3D Lennard- 
Jones gas from top to bottom. For large systems, p should 
scale as For small system sizes, the majority of 

particles is in the droplet p = 0{1) (horizontal dashed line). 
The Metropolis data are obtained in the droplet phase at 
the finite-size transition temperature as extrapolated from the 
linear fit in Fig. [5] 

the largest cluster p — or in other words the mass of the 
droplet divided by the total number of particles — that 
should scale as 

Figure [9] shows that this assumption seems to be valid 
for all three considered models. In each case, we show 
the expectation value (p)tc in the droplet phase at the 
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transition temperature. The horizontal dashed line at 
77 = 1 corresponds to all particles being in the droplet. 
For one, the results from our multicanonical simula¬ 
tion (red pluses) are measured as described in Sec. IIII El 
with jackknife errors. In addition, there are statistical 
averages from canonical Metropolis simulations (green 
crosses) prepared in the droplet phase at Tc{N) = Tq + 
as extrapolated in Sec. II V Al To be precise, 
we considered the following pairs (Tq, a): (0.623, —0.537) 
for the 3D lattice gas as well as (O.MO, —0.611) for the 3D 
Lennard-Jones gas. In this case, the errors are statisti¬ 
cal errors including the integrated autocorrelation times. 
The straight line is the expected scaling behavior shifted 
as a guide to the eye in order to be compared with the 
numerical data. 

The 2D lattice gas shows the expected scaling behav¬ 
ior quite early, already for roughly 100 particles. In the 
case of the 3D lattice gas the situation already changes 
and the expected scaling behavior only starts for quite 
large system sizes of the order of 2000 particles. This is 
consistent with the 3D Lennard-Jones gas, where, how¬ 
ever, the available system sizes are much smaller. Still 
the expected scaling may be anticipated. Including the 
Metropolis data, it shows that the leading-order extrap¬ 
olation of the available transition temperatures is not 
precise enough and seems to overestimate the transition 
temperature for larger systems, consistent with the ob¬ 
servations in Sec. IIV Al 

All cases show that for small system sizes, the majority 
of particles is included in the single macroscopic droplet. 
A corresponding constant behavior (7 ))to ~ 1 (the hori¬ 
zontal dashed line) would explain an intermediate scal¬ 
ing regime where the relevant length scale depends on the 
system volume. From Fig.|9l this may be at most guessed 
locally. However, regarding the measurement of the av¬ 
erage critical droplet we have to mention that the small 
systems including only a few particles also show narrow 
distributions. These are discretized by nature and only 
defined in the range [1,77]. Thus, the observable may 
not provide the required measurement sensitivity. Still, 
having a large fraction of particles in the droplet, espe¬ 
cially compared to the vanishing fraction in the limit of 
large systems, suffices for the argument that locally an 
effective intermediate scaling becomes visible. Moreover, 
it becomes clear that there are corrections to the size of 
the critical droplet for small and intermediate systems. 
This should also have an effect on the scaling of the tran¬ 
sition temperature and rounding. 

V. CONCLUSIONS 

We presented results on the canonical finite-size scal¬ 
ing of the droplet condensation-evaporation transition. 
To this end, we performed a comparative study in which 
we considered two- and three-dimensional lattice gas as 
well as three-dimensional off-lattice Lennard-Jones gas 
models in order to draw general conclusions. In all 
cases, we could verify the theoretical predictions about 


the finite-size corrections to the transition temperature 
oc and the finite-size rounding oc in 

the limit of large systems [ 1 , @ . In addition, we identified 
an intermediate regime with an effective different scaling 
behavior in all three cases, most apparent in the finite- 
size rounding. This regime is consistent with results from 
flexible polymers 0 , where the aggregate included most 
of the polymers. Our measurement of the average droplet 
size in equilibrium with vapor at coexistence showed that 
indeed for small system sizes a large fraction of the sys¬ 
tem is included in the single macroscopic droplet. The 
scaling of the droplet size, however, is nicely described 
by the theoretical predictions ia in the limit of large 
systems. 

Without knowledge of the larger system sizes, the in¬ 
termediate regime could give way to a wrong finite-size 
scaling behavior with qualitatively good fit results for 
both the transition temperature and rounding. However, 
we were able to show that there is a gradual crossover 
to the large-system scaling regime which can be approxi¬ 
mately described by including empirical higher-order cor¬ 
rections to the theoretical scaling predictions, namely 
polynomial orders of the critical droplet size. If only 
the leading-order corrections are considered, the round¬ 
ing is a good source getting a lower fit limit for a qual¬ 
itatively good fit. This was compared to analytic and 
low-temperature results for the lattice cases. It is thus 
also suitable in order to identify the finite-size scaling 
regime for phase transitions with mixed phases, like 
the condensation-evaporation transition but also more 
complex systems like polymer aggregation. Considering 
only the local, intermediate regime and applying stan¬ 
dard Hnite-size scaling approaches (with possibly wrong 
corrections) would then lead to wrong estimates of the 
infinite-system limit. A similar note of caution has been 
recently demonstrated in the finite-size scaling of self¬ 
avoiding walks on percolation clusters [s^j . 

An intuitive approach to the leading-order finite-size 
scaling corrections is the competition of volume {L‘^) 
and surface (L^“^) contributions. For a first-order phase 
transition, this gives rise to a finite-size correction of the 
order L~^, where L is the relevant length scale of the 
system 0, mu. We argue that the linear extension 
of the droplet R at coexistence plays this dominant role. 
The condensation-evaporation transition clearly connects 
a homogeneous and a mixed phase in the canonical en¬ 
semble. However, within the canonical approach, we may 
consider a (virtual) subsystem with the volume of the 
transition droplet. By translational invariance, this sub¬ 
system may be always constructed around the largest 
droplet. Above the transition temperature, this subsys¬ 
tem includes a homogeneous gas phase while at and below 
the transition it is filled by the largest droplet and hence 
shows a homogeneous liquid phase. Thus, this may be 
interpreted as a grand-canonical transition between ho¬ 
mogeneous phases in the virtual system spanned by the 
volume of the critical droplet. By construction, this vir¬ 
tual volume would have open boundary conditions yield¬ 
ing a finite-size shift of order R~^ and a finite-size round- 
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iiig of order R~'^. This picture is consistent with rigorous 
results for non-periodic first-order phase transitions [13 ■ 
However, this argument relies on the finite-size scaling 
of i? oc ^ which was shown to be already non¬ 

trivial. 

The present study shows that considering an orthogo¬ 
nal crossing of the phase boundary still yields the same 
finite-size corrections and serves as a complementary 
tool. On the example of the transition rounding, the 
fixed-density approach was shown to be closer to the ex¬ 
pected large-system scaling behavior already for smaller 
system sizes. In general, both directions have their ad¬ 
vantages and drawbacks both numerically and systemati¬ 
cally. This may be exploited for one’s benefit by choosing 
the suitable direction for the problem at hand, as was also 
recently demonstrated for the Blume-Capel model 
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